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ABSTRACT 

In  the  fall  of  1984  an  experiment  was  conducted  at  an  arctic  ice  camp  to  examine 
20-80  kHz  acoustic  reflections  from  an  isolated  ice  block.  The  results  are  compared  with 
a  simple  model  of  ice  block  reflections  that  has  been  used  in  conjunction  with  an  ice 
block  configuration  for  pressure  ridge  keels.  The  study  of  the  reflections  from  these  keels 
is  important  because  they  form  an  interfering  background  for  acoustic  equipment  operat¬ 
ing  under  the  ice.  A  cylindrical  block  of  ice  0.84  m  in  diameter  was  cut  from  the  flat  sur¬ 
face  of  a  floe  and  depressed  so  that,  for  a  transducer  placed  below,  the  reflection  from  the 
block  would  arrive  before  that  from  the  surrounding  ice.  A  transducer  15-30  m  below  the 
block  was  moved  horizontally  in  steps  to  measure  changes  with  aspect.  The  returns  from 
the  block  are  compared  with  those  predicted  by  the  model,  which  includes  a  loss  based 
on  the  bulk  impedance  properties  of  the  two  media.  When  near  field  effects  are  included, 
the  measured  returns  at  20  kHz  are  similar  to  those  predicted;  but  as  the  frequency  is 
increased  the  return  at  normal  incidence  is  lower  than  predicted  and  the  response  pattern 
is  broadened  and  smoothed.  The  complex  structure  of  the  growing  sea  ice,  producing 
both  volume  and  surface  acoustic  scattering,  is  the  most  likely  cause  of  the  discrepancy. 
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I.  INTRODUCTION 

In  the  arctic  the  ice  canopy  is  stressed  by  the  wind,  causing  differential  movements 
which  form  pressure  ridges  and  keels  of  all  shapes  and  sizes.  These  keels  are  a  disruptive 
influence  on  the  transmission  of  sound  beneath  the  ice  and  thus  have  a  large  effect  on  the 
performance  of  acoustic  equipment  in  this  medium. 

Measurements  of  ice-keel  reflections  in  1982  revealed  that  they  come  from  indivi¬ 
dual,  point-like  reflectors  which  tend  to  appear  and  disappear  with  small  changes  in 
aspect,  giving  a  sporadic  distribution  for  reflections  along  the  keel.1  Possible 
configurations  that  could  produce  such  reflections  include  blocks  with  a  face  toward  the 
source,  blocks  forming  a  comer  reflector,  and  blocks  that  are  tilted  so  as  to  reflect  to  the 
surface  and  back  to  the  receiver. 

In  the  fall  of  1984,  an  experiment  was  devised  to  examine  one  of  these 
configurations,  a  block  of  ice  with  a  face  at  near-normal  incidence.  For  this  experiment, 
a  cylindrical  block  was  cut  from  a  flat  section  of  a  floe  and  depressed  so  that,  for  a  trans¬ 
ducer  placed  below,  the  reflection  from  the  block’s  face  would  arrive  before  the 
reflection  from  the  under  surface  of  the  floe  and  thus  be  distinguishable.  As  the  echo 
amplitudes  were  recorded,  the  transducer  was  moved  successively  along  a  line  to  change 
the  aspect. 

The  immediate  purpose  of  this  report  is  to  compare  the  measured  falloff  of  target 
strength  at  incident  angles  near  normal  with  a  theoretical  approximation  presently  used 
for  ice  keel  modeling.2  This  theory  uses  the  Kirchhoff  approximation  to  combine  the 
change  in  impedance  at  the  water-ice  interface  and  the  interference  pattern  for  reflection 
from  a  circular  plate.  Near-field  effects  are  also  computed. 

The  results  show  that  at  20  kHz  the  measured  reflection  is  a  few  decibels  lower  than 
theory,  but  the  response  pattern  is  approximately  as  predicted.  At  higher  frequencies  the 
return  at  normal  incidence  is  considerably  lower  than  predicted,  and  the  measured 
response  pattern  drops  off  with  very  little  sign  of  the  nulls  and  peaks  that  are  present  in 
the  theoretical  diffraction  pattern.  Apparently  the  surface  and  volume  properties  of  the 
porous  ice  at  and  near  the  interface  provide  a  transition  layer  that  cannot  be  treated  as  a 
sharp  water-ice  boundary.  The  effect  of  the  gradual  transition  is  to  weaken  and  smooth 
the  response.  Further  investigation  is  needed  on  the  surface  and  volume  scattering  pro¬ 
perties  of  growing  sea  ice.  The  following  summary  of  the  freezing  process  may  be  help¬ 
ful. 


APL-UW  8506  1 


UNIVERSITY  OF  WASHINGTON  •  APPLIED  PHYSICS  LABORATORY 


Sea  ice  is  a  very  inhomogeneous  material.  Freezing  starts  at  the  sea  surface  with 
small  disk-shaped  particles  (platelets)  which  grow  horizontally  in  a  snowflake-like  pat¬ 
tern.  As  freezing  continues  and  the  platelets  begin  to  crowd  each  other,  they  eventually 
tip  to  a  vertical  position.  In  the  process,  brine  is  concentrated  between  them.  The  freez¬ 
ing  process  continues  with  the  downward  growth  of  the  ice  platelets.  Eventually  the 
space  between  the  platelets  freezes,  accompanied  by  displacement  and  some  entrapment 
of  the  brine  concentrate.  The  vertical  platelets  at  the  freezing  interface  form  a  "skeletal 
layer"  which  is  quite  fragile.  The  thickness  of  the  skeletal  layer  depends  on  the  tempera¬ 
ture  gradient  in  the  ice;  the  temperature  varies  from  about  -1.8°C  at  the  lower  surface  to 
the  air  temperature  at  the  upper  surface.  The  entrapped  brine  remains  in  solution,  in  tem¬ 
perature  equilibrium,  down  to  about  -23°C.  Below  -23°C,  most  of  the  brine  is  a  solid 
precipitate.  Above  -23°C,  for  a  given  bulk  salinity,  the  brine  volume  depends  on  the 
temperature.  Entrapped  air  and  gases,  displaced  from  seawater,  also  contribute  to  the 
porosity  of  ice  and  may  be  an  important  element  in  the  acoustic  response  of  sea  ice.  For 
the  ice  samples  studied  in  this  experiment,  the  close-off  of  pores  due  to  temperature 
probably  occurred  2  or  3  cm  from  the  platelet  tips,  giving  rise  to  an  "equivalent  acoustic 
roughness"  of  many  millimeters. 

The  result  of  the  experiment  will  be  helpful  in  determining  if  a  flat  ice  block  face  at 
normal  incidence  is  a  major  source  of  the  sporadic  reflections  from  an  ice  keel.  An 
identification  of  the  ice  configurations  in  a  keel  that  cause  acoustic  reflections  is  essential 
for  preparing  models  of  the  ice  for  acoustic  studies. 
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II.  EXPERIMENTAL  ARRANGEMENT 

The  experiment  was  c  onducted  on  an  ice  floe  in  the  Beaufort  Sea  on  3-5  November 
1984.  The  arrangement  for  the  experiment  is  shown  in  Fig.  1.  Polypropylene  lines  run¬ 
ning  through  oil-filled  pipes  inserted  through  the  ice  supported  the  transducer  from  three 
directions  120°  apart.  These  lines  were  tied  to  a  ring  above  the  surface,  and  the  trans¬ 
ducer  was  moved  by  moving  the  ring.  The  arrangement  on  the  surface  is  shown  in 
Fig.  2.  A  framework  near  the  submerged  ice  block  supported  a  beam  on  which  pegs 
were  spaced  every  15  cm.  All  four  sides  of  the  frame  were  notched  at  15  cm  intervals  so 
that  the  beam  containing  the  pegs  could  be  moved  sideways  or  rotated  perpendicularly. 
By  moving  the  ring  successively  along  the  line  of  pegs  above  the  ice,  the  transducer  was 
moved  along  a  corresponding  line  below  the  ice  with  an  exaggeration  that  was  easily  cal¬ 
culable.  The  calculation  was  facilitated  by  orienting  the  framework  so  that  its  centerline 
extended  through  one  of  the  pipes  holding  the  lines. 

Installation  began  with  divers  examining  the  under  side  of  the  floe  and  selecting  a 
relatively  flat  area.  A  4  in.  (10  cm)  diam.  hole  was  drilled  0.3  m  deep  in  the  0.46  m  thick 
ice  at  the  selected  location.  A  2  in.  (5  cm)  diam.  aluminum  pole  was  installed  as  verti¬ 
cally  as  possible  in  the  hole,  and  fresh  water  was  poured  in  to  help  freeze  the  pole  in 
place.  The  next  day,  a  0.84  m  diam.,  circular  block  was  cut  out  aro  ind  the  pole,  pushed 


Figure  1.  Arrangement  for  the  ice  block  experiment. 
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Figure  2.  The  framework  for  controlling  the  transducer  from  the  surface. 


downward  0.6  m,  and  guyed  into  position.  A  circular  block  was  used  so  that  its  orienta¬ 
tion  with  respect  to  the  frame  did  not  need  to  be  determined.  The  divers  then  helped 
install  the  lines  for  the  transducer.  Once  the  lines  were  available  at  the  surface,  the  trans¬ 
ducer  was  connected  and  lowered  to  the  desired  depth. 

Photographs  of  the  arrangement  as  installed  in  the  field  are  shown  in  Fig.  3.  Fuel  oil 
was  poured  into  each  of  the  three  plastic  pipes  through  which  the  lines  were  threaded  to 
keep  them  free.  The  electronics  cable  from  the  transducer  extended  100  m  to  a  heated 
building  which  housed  the  instrumentation.  Communication  between  the  field  crew  and 
the  instrument  room  was  by  radio. 
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|  Figure  3(b).  The  framework  placed  near  the  block.  A  2-in.  (5-cm)  diam.  aluminum 

tube  was  installed  vertically  in  the  ice  before  cutting  out  the  block. 

I 
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Figure  3(c). 

The  framework  placed  near  the 
block.  A  metal  ring,  placed  over 
one  of  the  pegs,  held  the  ends  of 
the  three  lines  at  the  desired 
location. 


Figure  3(d). 

Transducer  assembly  being 
lowered  through  the  diver 
access  hole.  A  lead  weight  kept 
the  transducer  vertical. 


6  APL-UW  8506 


UNIVERSITY  OF  WASHINGTON  •  APPLIED  PHYSICS  LABORATORY 


III.  THEORETICAL  TREATMENTS 

The  reflection  of  plane  waves  from  an  ice  block  face  has  been  predicted  by  Ellison2 
based  on  a  combination  of  two  separate  phenomena  which  are  incorporated  in  the 
Kirchhoff  approximation.3  One  is  the  reflection  loss  due  to  the  change  in  acoustic 
impedance  at  the  water-ice  boundary.  The  other  is  the  interference  pattern  produced  in 
the  reflected  wave  by  individual  contributions  from  separate  elements  of  the  finite  block 
face.  In  addition,  we  consider  near-field  effects  and  refer  to  ice  absorption  measurements. 


A.  Impedance  Change 

For  a  plane  wave  impinging  on  an  infinite  surface,  the  amplitude  of  the  reflection, 
Ar,  is  related4  to  the  amplitude  of  the  incident  plane  wave,  A,,  by  the  amplitude 
reflection  coefficient 


pc 


p,  ci 
pc 


+  B 


-l 


0) 


When  applied  to  the  water-ice  interface  p,  and  ci  are  the  density  and  sound  speed  for  ice, 
and  p  and  c  are  the  density  and  sound  speed  for  water; 


where  0,  is  the  incident  angle.  For  normal  incidence,  B  =  1. 

The  values  used  for  these  properties  are  shown  in  Table  I.  The  values  for  water 
were  measured  in  the  field,5  whereas  those  for  ice  were  estimated  from  measurements 
reported  by  others6  for  arctic  ice. 


Table  /.  Acoustic  properties  of  ice  and  water. 


Ice 

Water 

p,  =  0.92  g/cm3 

p  =  1.022  g/cm3 

c,  =  3200  m/s 

c  =  1433  m/s 

Impedance  Ratio  = 

PiC‘  =  2944  =  2.01 
pc  1465 
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Inserting  these  values  in  Eq.  1,  we  obtain  for  normal  incidence  (fi=l) 


*a  =  71  =  0-34  ,  (2) 

which  represents  a  decrease  in  intensity  of  9.5  dB  for  the  reflected  sound.  The  shear 
wave  should  also  be  considered,  but  its  effect  is  very  small  for  angles  of  0°-15°,  the  range 
of  angles  encountered  in  the  experiment. 


B.  Rigid  Plate  Reflections 

If  we  consider  the  block  face  as  a  fixed,  rigid  plate  (following  the  Ellison  approach), 
the  theoretical  target  strength  of  a  return  at  incident  angle  0  is  given  by7 


TS  =  20  log 


'iMP) 

P 

v 


cos  p 


(dB), 


(3) 


where 

P 

X 

k 

a 


(4n  a) 

X 


sin9 


wavelength 
2k  /  X 


radius  of  plate. 


For  this  equation  to  be  valid,  there  should  be  no  appreciable  near-field  effects  (i.e., 
R  >>a2/X,  where  R  =  range)  and  the  Kirchhoff  approximation  for  a  flat  plate  must  be 
valid  (i.e.,  ka  »  1).  For  our  measurements,  the  minimum  range  was  14.5  m  and  the 
maximum  a2lX  was  9.8,  indicating  some  near-field  effect.  Our  minimum  ka  was  37,  well 
above  1. 


C.  Combined  Effects 

If  we  combine  the  effects  of  the  impedance  change  and  the  finite  plane  pattern,  we 
obtain  the  following  equation  for  the  target  strength  of  the  block  face: 

K  a2  2  7 1  (P) 


TS  =  20  log 


R , 


P 


COS0 


(4) 


where  RA  =  (Q-1)/(Q+1);  Q  is  the  ratio  of  the  impedance  of  the  ice  to  that  of  the  water; 
i.e., 
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Equation  4  predicts  that  the  return  will  drop  rapidly  as  the  angle  departs  from  nor¬ 
mal  incidence,  but  in  addition  there  are  interference  peaks  and  valleys  whose  spacing 
depends  on  the  frequency.  Figure  4  shows,  as  examples,  the  returns  at  three  frequencies 
from  a  circular  plate  of  radius  0.42  m  as  computed  from  Eq.  4. 
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Figure  4. 

Interference  patterns  pre¬ 
dicted  for  reflection  from  a 
rigid,  circular  plate  com¬ 
bined  with  the  loss  (95  dB) 
due  to  the  impedance  change. 


D.  Near-Field  Effects 

For  short  ranges,  the  possibility  of  near-field  effects  should  be  investigated.  The 
possibility  can  be  assessed  by  comparing  the  phase  of  contributions  from  the  edge  of  the 
block  face  with  that  of  those  from  the  center,  at  normal  incidence.  The  phase  difference 
P  is  calculated  by 
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P  =2 


^  R2  +  a2-R]  ^ - 

J  X 


(radians) , 


where  a  is  the  radius  of  the  face.  Table  II  shows  the  phase  differences  calculated  for  the 
frequencies  and  ranges  used  in  the  experiment.  For  the  shorter  range,  the  phase  shift  is 
considerable  at  40  kHz  and  higher. 


Table  11.  Phase  difference  at  edge  of  face. 


Frequency 

(kHz) 


Phase  Difference  (deg) 

R  =  14.5  m  R  =  32.9  m 


61 

27 

122 

54 

47 

143 

63 

244 

107 

(radius  a  =  0.42  m) 

ra-i 


To  determine  the  effect  of  the  phase  differences,  we  numerically  integrate  outward 
over  concentric,  circular  rings,  adding  the  in-phase  component  of  each  area.  The  sum 
divided  by  the  area  is  designated  / .  Calculated  values  of  /  are  given  in  Table  III  along 
with  the  calculated  reduction  in  target  strength,  -20  log/.  Table  HI  shows  that  the 
short-range  measurements  require  an  appreciable  near-field  correction. 


Table  III.  Near-field  effect. 


Target  Strength  Reduction  (dB) 

Frequency  Area  Factor/  (Runs  5-7)  (Runs  1-3) 

(kHz)  R  =  14.5  m  R  =  32.9  m  R  =  14.5  m  R  =  32.9  m 


0.95 

0.99 

0.4 

0.83 

0.96 

1.6 

47 

0.77 

0.95 

2.3 

0.4 

80 

0.86 

7.7 

1.3 
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The  angular  dependence  can  be  computed  by  numerically  integrating  over  ring  seg¬ 
ments.  The  results  for  the  shortest  range  used  in  the  experiment,  15  m,  are  shown  in 
Fig.  5.  The  main  lobe  is  reduced,  and  the  nulls  at  the  smaller  angles  tend  to  be  bridged 
across. 


Figure  5. 

Computed  target  strength 
with  corrections  for  near¬ 
field  effects  at  a  range  of 
15  m.  Compare  with  Fig¬ 
ure  4. 
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E.  Sound  Absorption  in  Ice 

The  simple  model  proposed  by  Ellison  does  not  account  for  any  contribution  to  the 
observed  target  strength  from  volume  scattering  or  from  rear  surface  reflection.  In  our 
analysis,  we  investigated  a  second  reflection,  apparently  from  the  upper  surface  of  the 
block,  which  was  originally  an  air-ice  interface.  Evaluating  the  returns  from  the  upper 
face  of  the  block  requires  knowledge  of  the  absorption  coefficient  in  the  ice.  With  no 
theory  available,  we  used  a  recent  summary  of  absorption  measurements  by  McCammon 
and  McDaniel6  which  indicates  that  the  absorption  of  compressional  sound  waves  of  fre¬ 
quency  /  (kHz)  in  ice  of  temperature  T  (°C)  is 

a  =  0.06/  (-6/T)2'3  (dB/m)  (6) 

for  T  limited  to  -2°  to  -20°C. 

The  air  temperature  was  about  -20°C  during  the  days  prior  to  cutting  the  block.5 
The  water  temperature  was  -1J°C.  After  the  block  was  submerged  for  4  days,  its  aver¬ 
age  temperature  was  somewhere  between  the  initial  average  of  -1 1°C  (-20°  at  the  top  and 
-1.5°C  at  the  bottom)  and  the  water  temperature  of  -1.5°C.  Table  IV  shows  the  absorp¬ 
tion  for  the  temperatures  and  frequencies  of  interest  as  calculated  using  Eq.  6. 


Table  IV.  Absorption  of  sound  in  ice 6  (decibels/meter)  for  temperatures  between  -20°C 
and  -2°C,  as  calculated  using  the  equation  a  =  0.06/  (-6/T  )2/3  dB/m. 


T(°C) 

20 

Frequency  (kHz) 

40  47 

80 

-1.5° 

3.0 

6.1 

7.1 

12.1 

-5 

1.4 

2.7 

3.2 

5.4 

-10 

0.9 

1.7 

2.0 

3.4 

-15 

0.7 

1.3 

1.5 

2.6 

-20 

0.5 

1.1 

1.3 

2.2 

a Extrapolated  slightly  beyond  limit  of  equation. 

NOTE:  The  two-way  path  through  the  ice  block  was  0.92  m. 
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IV.  MEASUREMENTS 

To  position  the  transducer  relative  to  the  ice  block,  a  gridwork  was  established  on 
the  surface  with  position  0,0  in  the  southwest  comer.  The  positions,  which  were  at  15  cm 
intervals,  were  labeled  according  to  their  distance  north  or  east  of  this  point.  Each  coor¬ 
dinate  had  32  positions.  A  wooden  beam  with  vertical  pegs  spaced  every  15  cm  was 
oriented  north-south  or  east- west  in  a  square  framework  notched  at  15  cm  intervals.  In 
this  way,  any  line  of  peg  positions  could  be  provided.  A  set  of  measurements  labeled 
line  N20,  for  example,  would  be  obtained  with  the  beam  oriented  east-west  in  the  20th 
position  north  of  0,0.  The  ring  connected  to  the  lines  holding  the  transducer  would  then 
be  moved  along  the  line  from  one  peg  to  another. 

For  the  experiment,  short  pulses  were  transmitted  at  either  20,  40,  47,  or  80  kHz  and 
received  on  the  same  transducer.  The  returns  were  observed  for  about  10  s  on  a  Nicolet 
digital  oscilloscope,  and  the  average  amplitude  of  the  reflections  was  estimated.  After 
several  trials  to  locate  the  position  producing  the  maximum  return,  the  ring  was  moved 
from  peg  to  peg  along  two  perpendicular  lines  intersecting  that  position,  and  the  average 
amplitude  at  each  location  was  recorded  by  hand.  In  addition,  one  set  of  returns  at 
80  kHz  was  digitized  on  the  oscilloscope,  and  the  results  were  recorded  on  a  diskette. 
This  made  possible  later  examination  of  the  echo  in  detail. 

The  center  of  the  framework  was  at  grid  position  N  16-El 6,  but  this  was  not  neces¬ 
sarily  the  location  that  placed  the  transducer  directly  beneath  the  block.  On  the  first  day 
the  center  of  the  block  appeared  to  be  at  N  18-El  8.  After  the  lines  were  shortened  for  the 
measurements  on  the  second  day,  the  center  was  near  N20-E18. 

A.  Data  Summary 

The  first  three  "runs"  were  made  on  3  November  with  the  transducer  at  a  depth  of 
33.9  m.  At  that  depth,  the  transducer  moved  horizontally  41%  more  than  the  ring.  On 
5  November  the  lines  were  shortened,  raising  the  transducer  to  a  depth  of  15.5  m,  and 
three  more  runs  were  made.  At  this  depth,  the  transducer  moved  only  about  10%  more 
than  the  ring,  and  the  returns  were  stronger.  The  runs  are  summarized  in  Table  V.  The 
pulse  length  was  1.0  ms  for  Runs  1  and  2,  and  0.5  ms  for  the  remaining  runs. 

B.  Two  Echoes  in  the  Return 

The  envelopes  of  the  returns  on  Run  6,  which  used  an  80  kHz,  0.5-ms-long  cw 
pulse,  were  digitized  on  the  Nicolet  oscilloscope,  preserved  on  diskette,  and  examined 
later  in  detail.  A  sample  return  is  shown  in  Fig.  6.  Note  the  two  overlapping  echoes.  The 
second  echo  arrived  only  0.3  ms  after  the  first  and  well  before  the  echo  from  the  under 
side  of  the  adjacent  floe,  the  start  of  which  is  shown  at  the  right  in  the  figure.  The 
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Table  V.  Measurements  of  reflections  from  an  ice  block. 


Run 

No. 

Date 

Transducer 

DepthJ 

(m) 

Frequency 

(kHz) 

Transmitter 

Voltage* 

(V) 

Total 

Gain 

(dB) 

Type  of 
Transducer 

Receiving 

Sensitivity 

(dB)c 

Transmitting 

Response 

(dB)d 

1 

3  Nov  84 

33.9 

47 

200 

79 

2"  piston 

-195 

156 

2 

47 

200 

79 

-195 

156 

3 

5  Nov  84 

80 

620 

80 

ITC  1042 

-207 

150 

5 

15.5 

40 

660 

82 

-208 

136 

6 

80 

530 

82 

-207 

150 

7 

20 

520 

94 

-207 

123 

“The  lower  face  of  the  block  was  at  a  depth  of  1  m. 
*  Peak-to-peak 

c  Decibels  are  referenced  to  1  V/|iPa. 

^  Decibels  are  referenced  to  1  jiPa/V  at  1  m. 


Figure  6. 

Return  from  the  ice  block  on  Run  6  at 
80  kHz,  N20-EI6,  0.5  ms  cw  pulse. 


dimensions  of  the  arrangement  (Fig.  7)  and  the  expected  sound  speed  in  the  ice  enable  us 
to  calculate  that  the  first  echo  is  from  the  lower  face  of  the  block,  the  second  echo  is  from 
the  upper  face,  and  the  following  echoes  are  from  the  adjacent  ice.  Figure  8  compares 
the  levels  of  the  first,  second,  and  overlapping  returns  for  line  N20  of  Run  6.  The  overlap 
should  represent  the  vector  sum  of  the  echoes  from  the  top  and  bottom  faces  of  the  block. 
The  amplitude  of  the  overlap  was  near  the  algebraic  sum  for  the  two  echoes,  indicating 
that  the  amplitudes  added  coherently  for  the  particular  thickness  of  the  ice  block  used  in 
the  experiment.  The  phase  difference  would  include  a  1 80°  change  at  the  upper  face. 
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Figure  7.  Spacing  of  the  ice  block  for  the  experiments  on  3-5  November  1984. 
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Figure  8.  Comparison  of  first  echo,  second  echo,  and  overlapping  of  the  two  in  the 
return  at  80  kHz  on  Run  6;  line  N20. 
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The  presence  of  echoes  from  the  upper  face  at  the  other  frequencies  is  uncertain. 
Two  sample  oscillograms  of  the  returns  on  Run  2,  which  used  a  47  kHz,  1  ms  long  pulse, 
were  saved,  and  they  show  a  gradual  rise  in  amplitude  rather  than  distinct  levels.  Because 
of  experiment  time  constraints,  no  oscillograms  were  saved  of  the  measurements  at  20 
and  40  kHz,  which  used  a  0.5  ms  pulse.  Furthermore,  an  additional  reflection  from  the 
upper  face  of  the  ice  block  was  not  readily  apparent  in  any  of  the  echoes  examined  in 
detail  during  the  experiment.  Readings  of  maximum  signal  level  were  made  in  real  time 
and  logged  manually;  consequently,  the  measurements  could  represent  the  maximum  of 
the  first,  second,  or  overlapping  echo.  For  the  detailed  observations  at  80  kHz,  which  did 
show  two  echoes,  the  difference  between  the  amplitude  of  the  first  echo  and  the  max¬ 
imum  amplitude  was  small  (Fig.  8).  The  smaller  absorption  at  the  lower  frequencies 
would  tend  to  enlarge  the  returns  from  the  upper  surface,  but  with  no  multiple  levels 
observed  we  are  uncertain  as  to  the  effect  of  the  second  surface  on  our  readings. 

C.  Returns  at  Normal  Incidence 

The  80  kHz  echoes  were  examined  in  the  returns  at  incident  angles  near  zero.  The 
amplitudes  of  returns  at  these  angles,  which  were  assumed  to  represent  the  lower  face, 
are  tabulated  in  Table  VI  along  with  the  target  strength  calculated  from  the  transducer 
calibrations  and  the  two-way  transmission  losses.  Also  shown  are  the  reflection 
coefficients,  RA ,  computed  by  solving  for  RA  in  Eq.  4  and  making  the  near-field  correc¬ 
tion  given  in  Table  III. 

Equation  2  predicts  an  RA  of  0.34  for  the  density  and  sound  speed  values  involved 
in  the  experiment.  The  value  shown  in  Table  VI  for  20  kHz  (0.24)  is  a  little  below  0.34, 
while  the  values  shown  for  the  other  frequencies  are  far  below,  indicating  that  Eq.  4  is 
unacceptable  at  these  frequencies.  An  additional  correction  involving  the  scattering 
processes  at  the  lower  surface  is  apparently  needed,  the  more  so  the  smaller  the 
wavelength. 

As  mentioned  earlier,  the  reflection  measurements  on  Run  6  at  80  kHz  gave  the 
amplitude  of  the  return  from  the  upper  face.  The  reflection  coefficient  for  the  upper  face 
can  be  computed  by  comparing  the  reflections  from  the  upper  and  lower  surfaces,  if  some 
rather  bold  assumptions  are  made.  We  assume  that  the  reflection  coefficient  at  the  lower 
face  is  given  by  the  impedance  ratio  (Eq.  2,  RA  =  0.34)  and  that  both  returns  are  affected 
the  same  by  the  shape  of  the  ice  face.  Reflection  from  the  lower  face  changes  the  ampli¬ 
tude  of  the  incident  sound  by  the  factor  RA .  For  sound  entering  the  ice,  the  amplitude  is 
changed  by  1+/?A ;  at  the  reflection  off  the  upper  face  it  is  changed  by  an  unknown  factor, 
RB ,  and  as  it  passes  through  the  lower  face  it  is  changed  by  1  -RA .  In  addition,  the  sound 
suffers  absorption  losses  both  ways  through  the  ice.  These  changes  are  shown  in  the  fol¬ 
lowing  equation  for  the  difference  in  the  two  returns: 

C TS  )B  -  C TS  )A  =  20  loga-i-^ )  +  20  log(l-/?* )  +  20  log RB  -  20  log/?A  -2a  r  ,  (7) 
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Table  VI.  The  echo  amplitudes  used  to  calculate  target  strength,  and  the  amplitude 
reflection  coefficients  RA  computed  for  the  lower  face ,  for  all  runs. 

Reflection  Target 

Frequency  Transducer  Amplitude  Strength  RA  from 

Run  (kHz)  Position  (V,  peak)  (dB)  Average  Eq.  4 
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where  t  is  the  thickness  of  the  ice  block.  For  a  given  RA  and  a,  the  equation  can  be 
solved  for  RB ,  the  reflection  coefficient  for  the  upper  face. 

The  absoiption  in  the  ice  varies  with  temperature  as  shown  in  Table  IV.  The  tem¬ 
perature  of  the  ice  was  initially  -20°C  at  the  upper  surface  and  -1.5°C  at  the  lower.  Run  6 
was  made  after  the  ice  had  been  submerged  for  2  days.  The  average  temperature  was 
somewhere  between  the  initial  average,  -11°C,  and  an  ultimate  equilibrium  temperature 
of  -1.5°C.  In  calculating  RB,  the  reflection  coefficient  for  the  upper  surface,  we  used  both 
-6  and  -2°C  to  show  how  the  calculated  RB  depends  on  the  assumed  temperature. 

The  results  are  plotted  in  Fig.  9  for  two  lines  of  Run  6,  N20  and  E18.  The  reflection 
coefficient  depends  strongly  on  the  value  used  for  the  sound  absorption  in  the  ice,  vary¬ 
ing  by  a  factor  of  2  for  the  given  variation  in  ice  temperature.  The  coefficients  for  -6°C, 
the  more  likely  average  temperature,  are  near  the  value  (0.34)  calculated  for  the 
impedance  change  from  ice  to  water,  indicating  that  the  upper  surface  (originally  snow 
covered)  represents  the  same  change  in  density  and  sound  speed  as  the  lower  face. 

When  the  reflections  from  the  upper  and  lower  surfaces  overlap,  the  result  appears 
to  be  a  coherent  addition  of  the  two  reflections  with  a  phase  difference  of  about  45°  as 
shown  in  the  upper  graph  of  Fig.  9.  The  amplitude  of  the  overlap  is  definitely  greater  than 
for  an  incoherent  addition;  coherence  seems  to  prevail  even  after  the  sound  has  passed 
through  the  block  and  reflected  from  the  upper  surface. 

At  the  other  frequencies,  the  contribution  from  the  upper  face  was  not  detected  dur¬ 
ing  the  experiment,  and  no  records  are  available  for  further  analysis. 


Figure  9. 

Upper  surface  reflection  coefficients 
calculated  for  Run  6  at  80  kHz  ( bot¬ 
tom )  along  with  calculated  relative 
phase  of  the  two  reflections  (top). 
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D.  Aspect  Dependence 

To  investigate  the  effect  of  changes  in  aspect  angle,  the  transducer  was  moved  suc¬ 
cessively  along  a  line  beneath  the  block.  At  each  position,  the  location  of  the  transducer 
and  the  aspect  angle  of  the  sound  incident  on  the  block’s  face  were  calculated  from  the 
geometry  of  the  supporting  lines.  The  amplitude  of  the  return  was  then  converted  to  tar¬ 
get  strength  and  plotted  versus  the  calculated  aspect  angle.  The  results  are  shown  in 
Fig.  10.  The  characters  following  the  run  number  indicate  the  line  along  which  the  ring 
controlling  the  transducer’s  position  was  moved. 

The  data  points  in  Fig.  10  are  compared  with  the  theoretical  pattern  predicted  by 
Ellison,2  who  uses  the  Kirchhoff  approximation  in  combining  impedance  change  and 
rigid  plate  reflection.  Near-field  effects  are  included  in  the  figure  for  the  short  range  runs 
(5,  6,  and  7).  The  high  peak  predicted  near  0°  was  never  observed,  but  the  dropoff  with 
angle  beyond  ±2°  is  about  as  predicted.  There  is  an  occasional  indication  of  the  predicted 
side  lobes,  but  the  minima  are  spaced  farther  apart.  Some  lines  did  not  pass  directly 
below  the  ice  block;  for  plots  of  these,  the  theoretical  return  was  calculated  for  an  offset 
from  center.  The  predicted  sharp  minima  were  detected  so  rarely  that  at  this  time  they 
will  be  ignored. 


Figure  10. 

Response  patterns  for  a 
0.84  m  diam.  ice  block.  The 
lines  are  the  theoretical  pat¬ 
tern  for  a  rigid  plate,  com¬ 
bined  with  impedance 
change.  Corrections  for  near 
field  have  been  made  for  the 
short  range  measurements. 
(Continued  on  the  following 
pages.) 
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Figure  10,  corn. 


INCIDENT  ANGLE  (deg) 
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Figure  10,  cont. 
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Figure  10,  cont. 


INCIDENT  ANGLE  (deg) 


To  test  the  form  of  the  aspect  dependence,  the  plots  of  target  strength  versus 
incident  angle  in  Fig.  10  were  folded  over  to  combine  plus  and  minus  incident  angles, 
and  the  data  were  fitted,  using  a  least-squares  method,  by  the  line 

TS=A-BQ.  (8) 

As  shown  in  Figs.  1  la-c,  the  lines  appear  to  represent  the  data.  Considering  the  scatter  in 
the  data,  a  more  complex  dependence  on  angle  is  not  warranted. 

The  constants  A  and  B  for  each  line  are  listed  in  Table  VII  along  with  the  average 
for  each  set.  The  averages  for  A  have  been  adjusted  for  the  near-field  correction  shown  in 
Table  III.  Note  that  the  80  kHz  data  for  the  two  transducer  depths  are  brought  into  agree¬ 
ment  by  the  near-field  correction. 
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Table  VII.  Constants  for  the  line  TS  =A  -  B  9. 


Frequency 

(kHz) 

Run 

Corrected6 

A  A 

Average 

Corrected 

A 

B 

Average 

B 

20 

7N20 

5.7 

6.1 

2.88 

7E18 

5.5 

5.9 

2.90 

7E15 

3.3 

3.7 

5.2 

3.05 

2.9 

40 

5N20 

-5.7 

-4.1 

2.25 

5E18 

-8.8 

-7.2 

-5.7 

1.44 

1.9 

47 

1N18 

-12.3 

-11.9 

3.10 

1E18 

-10.7 

-10.3 

3.15 

2N16 

-18.5 

-18.1 

3.29 

2N18 

-12.3 

-11.9 

4.23 

2N19 

-13.8 

-13.4 

3.66 

2N20 

-16.5 

-16.1 

3.16 

2E16 

-16.3 

-15.9 

4.05 

2E17 

-13.4 

-13.0 

4.41 

2E18 

-12.4 

-12.0 

4.43 

2E20 

-20.8 

-20.4 

-14.3 

2.15 

3.6 

80 

3N19 

-10.1 

-8.8 

2.41 

3N20 

-6.3 

-5.0 

2.90 

3N21 

-7.2 

-5.9 

3.05 

3E15 

-12.8 

-11.5 

2.69 

3E17 

-11.5 

-10.2 

3.06 

3E18 

-8.5 

-7.2 

2.64 

3E19 

-8.5 

-7.2 

-8.1 

3.10 

2.8 

6N20 

-14.2 

-6.5 

2.31 

6N20a 

-13.7 

-6.0 

2.37 

6E18 

-16.4 

-8.7 

1.83 

6E18a 

-17.4 

-9.7 

-7.7 

1.94 

2.1 

a  Using  first  echo  only. 

b Corrected  for  near-field  effect,  using  Table  III. 
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The  constant  A  is  plotted  against  frequency  in  Fig.  12  and  compared  with  the 
theoretical  prediction  based  on  the  impedance  change  and  rigid  plate  scattering.  The 
discrepancy  is  fairly  small  at  20  kHz,  but  exceeds  20  dB  at  the  higher  frequencies.  At 
these  frequencies,  it  is  also  more  or  less  constant.  Such  behavior  would  be  expected  if  the 
effective  roughness  of  the  under  surface  was  enough  to  cause  phase  changes  in  the  return 
greater  than  2n,  since  this  effect  should  reach  a  saturation  point  when  the  rms  roughness 
a  is  about  X/4.  At  47  kHz,  this  saturation  point  would  occur  at  o  =  0.8  cm,  which  is  in 
reasonable  agreement  with  the  size  of  the  skeletal  layer  observed  on  the  bottom  of  other 
blocks  cut  from  the  same  area. 


Figure  12.  Measured  target  strengths  at  normal  incidence,  i.e.,  values  of  A  in 
Table  VII,  corrected  for  near-field  effects. 


Although  the  dropoff  in  the  return  with  increasing  frequency,  and  the  eventual 
saturation,  appears  to  be  an  equivalent  roughness  effect,  we  do  not  have  enough  data  to 
suggest  an  equation  for  the  constant  A .  At  20  kHz  the  effect  is  small;  at  40-80  kHz  it 
may  be  20-30  dB. 
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A  plot  of  constant  B  versus  frequency  (Fig.  1 3)  for  the  five  sets  of  data  shows  con¬ 
siderable  variation  but  no  frequency  dependence. 


FREQUENCY  (kHz) 


Figure  13.  The  constant  B  in  Eq.  8  as  determined  from  the  measurements  listed  in 
Table  VII. 
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V.  SUMMARY 

The  results  of  the  1984  ice  block  experiment  can  be  summarized  as  follows: 

1.  The  returns  from  an  ice  face  at  near  normal  incidence  are  smaller  than  the 
theoretical  predictions  calculated  by  combining  the  rigid  plate  pattern  with  the 
impedance  change  effect  as  in  the  Kirchhoff  approximation.  The  difference  varies  from 
3  dB  at  20  kHz  to  as  much  as  30  dB  at  47  kHz  and  above,  and  may  be  due  to  the  inade¬ 
quacy  of  the  theory  or  to  the  structure  at  the  ice- water  boundary. 

2.  The  acoustic  response  of  the  0.84  m  diam  ice  block  dropped  off  about  10  dB  at 
4°  from  normal  incidence  for  all  frequencies  (20-80  kHz).  The  -3  dB  point  occurred  at  a 
beam  width  of  2.3°.  There  was  some  variation  in  the  beam  width  with  frequency  but  no 
noticeable  t  end.  Some  sharp  minima  were  observed  in  the  acoustic  response,  but  they 
did  not  agree  with  the  side  lobes  predicted  for  a  rigid  plate. 

3.  The  internal  reflection  from  the  submerged  upper  face  of  the  block  was 
observed  at  80  kHz,  indicating  that  the  reflections  from  the  back  side  of  an  ice  block 
deserve  further  investigation. 
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VI.  DISCUSSION 

In  earlier  work,1,8-11  a  detailed  examination  of  some  of  the  reflections  from  ice 
keels  showed  them  to  come  from  individual,  separated  reflectors  with  a  narrow  response 
pattern.  A  study  of  ice  block  configurations  that  would  produce  such  reflections  is  neces¬ 
sary  to  understand  ice-keel  reflections. 

The  1984  measurements  reported  here  were  made  to  examine  a  configuration  that 
seems  a  likely  source  of  such  reflections,  an  ice  block  face  at  normal  incidence.  These 
may  be  the  first  measurements  with  an  isolated  ice  block.  Good  control  of  the  incident 
angle  was  important  to  the  experiment.  An  examination  of  the  bottom  of  the  block  by 
divers  was  desirable  but  prohibitive  because  of  the  possible  contamination  of  the  under 
surface  with  air  bubbles.  Additional  measurements  would  have  been  prudent  but  the  per¬ 
sonnel  and  facilities  were  required  for  other  projects. 

The  measurements  indicate  that  this  configuration  could  be  a  major  source  of  ice 
keel  reflections  observed  in  the  field,  for  the  following  reasons:  (1)  the  measured  target 
strength  of  the  return  (5  dB)  and  size  of  the  block  (0.84  m)  correspond  to  those  for  a  keel, 
(2)  the  response  pattern  was  narrow,  as  observed  in  studies  of  returns  from  ice  keels,  (3) 
rough  estimates  indicate  that  there  may  be  enough  block  faces  present  at  near-normal 
incidence  in  a  keel  to  produce  a  good  portion  of  the  observed  reflectors. 

By  itself,  the  Kirchhoff  approximation,  which  allows  combining  the  effect  of  the 
impedance  change  on  specular  reflections  of  plane  waves  with  the  interference  pattern  for 
a  rigid  plate  at  a  given  incident  angle,  appears  to  be  inadequate  for  predicting  the  return 
from  in-situ  ice  surfaces  at  near-normal  incidence.  The  effect  of  the  porous  ice  structure 
at  the  lower  boundary,  which  becomes  more  important  at  the  shorter  wavelengths,  must 
be  included.  For  the  shorter  wavelengths,  this  porous  layer  gives  rise  to  a  gradient  in  the 
characteristic  impedance.  At  longer  wavelengths,  the  effective  impedance  mismatch 
between  ice  and  water  will  be  more  abrupt,  and  thus  the  interface  is  more  reasonably 
modeled  by  considering  the  Kdk  properties  of  the  two  media. 

The  peak  in  the  response  pattern  at  normal  incidence  is  a  demonstration  of 
Huygens’  principle  which  gives  the  return  in  the  far  field  as  the  sum  of  the  in-phase  con¬ 
tributions  from  all  area  increments  of  the  surface.  Any  surface  irregularity,  if  at  an 
appreciable  fraction  of  a  wavelength,  will  greatly  reduce  this  central  peak,  as  was 
observed  in  the  measurements. 

The  discovery  that,  at  80  kHz,  the  back  side  of  the  block  contributes  to  the  reflection 
means  that  further  measurements  must  be  made  to  investigate  interior  reflections  and  the 
absorption  in  the  ice.  The  effect  of  the  skeletal  layer  on  the  lower  surface  of  the  ice  and 
of  the  submerged  snow-covered  surface  above  must  also  be  investigated. 
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Some  sharp  minima  were  detected  in  the  return  at  angles  off  normal  incidence,  but 
the  spacing  did  not  correspond  to  that  predicted  for  an  idealized  flat,  uniform  surface.  It 
would  be  helpful  in  the  future  to  include  a  test  on  a  block  with  the  bottom  layer  cut  off; 
the  pattern  of  the  return  from  a  flat  solid  surface  could  then  be  measured  in  detail.  This 
would  eliminate  the  skeletal  layer  roughness  and  give  a  test  of  the  use  of  the  Kirchhoff 
approximation  for  an  ice  face. 

Additional  measurements  are  planned  for  1986,  using  a  block  cut  from  thicker  ice. 
After  the  first  measurement,  the  block  will  be  removed,  the  bottom  layer  cut  off,  and  the 
block  reset  for  another  measurement.  This  will  show  the  effect  of  the  skeletal  layer  and 
give  a  measure  of  the  absorption  in  the  ice,  if  the  upper  surface  is  detectable.  As  a  refer¬ 
ence  for  calibration  of  the  system  we  will  also  replace  the  ice  block  with  a  cylindrical 
block  of  air  to  give  a  flat  surface  with  a  high  impedance  change.  Several  diameters  will 
be  needed  as  a  check  on  the  area  dependence  of  the  theory.  The  acoustic  measurements 
will  be  supplemented  by  detailed  studies  of  the  density,  salinity,  temperature,  and  void 
ratio  of  the  ice,  especially  near  the  newly  frozen  surface.  The  lower  surface  will  be 
characterized  as  to  roughness  during  several  stages  of  the  freezing  process. 
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